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ABSTRACT 

The stages before the formation of stars in molecular clouds are poorly understood. 
Insights can be gained by studying the properties of quiescent clouds, such as their 
magnetic field structure. The plane-of-the-sky orientation of the field can be traced 
by polarized starlight. We present the first extended, wide-held (~10 deg 2 ) map of 
the Polaris Flare cloud in dust-absorption induced optical polarization of background 
stars, using the RoboPol polarimeter at the Skinakas Observatory. This is the first 
application of the wide-held imaging capabilities of RoboPol. The data were taken in 
the R-band and analysed with the automated reduction pipeline of the instrument. 
We present in detail optimizations in the reduction pipeline specihc to wide-held ob¬ 
servations. Our analysis resulted in reliable measurements of 648 stars with median 
fractional linear polarization 1.3%. The projected magnetic held shows a large scale 
ordered pattern. At high longitudes it appears to align with faint striations seen in 
dust emission, while in the central 4-5 deg 2 it shows an eddy-like feature. The overall 
polarization pattern we obtain is in good agreement with large scale measurements by 
Planck of the dust emission polarization in the same area of the sky. 

Key words: stars: formation - ISM: clouds - ISM: individual objects (Polaris Flare) 
- polarization. 


1 INTRODUCTION 

Molecular clouds in their vast complexity hold the key to un¬ 
derstanding the early stages of the star formation process. 
Magnetic fields and turbulence are the two main mecha¬ 
nisms that dictate the structural, dynamical and evolution- 
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ary properties of these clouds, through their competition 
against gravity. Their role in the onset of star formation can 
be studied best in quiescent non-star-forming regions, where 
stellar feedback is not present. One such region is the Polaris 
Flare, a translucent molecular cloud at a d istance between 
130 pc and 2 40 pc dHeithausen et al~lfl9 93l. discovered by 
iHeiles I J1984T) . It is believed to be in the early stages of its 
formation, since it does not ex hibit any signs of active star 
formation (I Andre et al. lljoTol ). CO observations have pro¬ 
vided invaluable information on th e turbulence signatures 
in the densest parts of the cloud (IFalgarone et al.l Il998l : 
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iHilv-Blant fe Falgarone 1120091 ). Recently, the Herschel space 
telescope mapped over 15 deg 2 of the cloud in dust emission 
dAndri et, al. Il201fil : iMiville-Deschenes et al. 11 2()1 f il l . 

The structure of the magnetic field of a cloud, as pro¬ 
jected on the plane of the sky, can be probed by observing 
polarized radiation. The polarization of starlight transmit¬ 
ted through a cloud is believed to be caused by dichroic 
extinction due to aspherical dust grains that are partially 
aligned with the magnetic field. This alignment causes the 
polarization direction of the light of background stars to 
trace the magnetic field direction of the cloud as projected 
on the plane of the sky. The same alignment process causes 
the thermal emission of these dust grains to be polarized in 
the direction perpendicular to the magnetic field. 

Information on the plane-of-the-sky magnetic field of 
the Polaris Flare has been pr ovided for the first time throu gh 
polar iz ed dust emission ([P lanck c ollaboration Int. XIX I 
120141 : IPlanck collaboration Int. XX l201ll ). These data, 
however, are limited by the instrumental resolution and con¬ 
fusion along the line of sight. A mapping of the region in po¬ 
larized starlight, which is complementary to the dust emis¬ 
sion but suffers from different limitations, is necessary to 
resolve these issues. 

We obtained optical polarization measurements of stars 
projected on 10 deg 2 of the Polaris Flare region with 
RoboPol. The RoboPol instrument is a 4-channel optical 
polarimeter with no moving parts, other than a filter wheel 
(Ramaprakash et al., in prep). It can measure both linear 
fractional Stokes parameters q = Q/I and u = U/I simulta¬ 
neously, thus avoiding errors caused by the imperfect align¬ 
ment of rotating optical elements and sky changes between 
measurements (polarization, seeing conditions). 

Each star in the field of view creates four images (spots) 
on the CCD displaced symmetrically in the horizontal and 
vertical directions. A mask supported by four legs is posi¬ 
tioned at the centre of the field of view. This allows targets 
that are centred on the mask to be measured with four times 
lower sky noise than the rest of the field. A typical image 
seen with RoboPol is shown in Fig. [T] The instrument has a 
13 r x 13 / field of view, enabling the rapid polarimetric map¬ 
ping of large areas of the sky. RoboPol is equipped with stan¬ 
dard Johnson-Cousins R- and /-band filters and is mounted 
on the 1.3-m, i/7.7 Ritchey-Cretien telescope at Skinakas 
Observatory in Crete, Greece. It has been operating since 
May 2013. 

RoboPol has been monitoring the optical linear polar¬ 
ization of a large s ample of gamma-ray bright blazars for 
the past two years (|Pavlidou et al.l [2014i j. In addition, the 
instrument is b eing used for lon g-term monitoring of Be 
X-ray binaries (IReig et al. 1120140 . Observations of optical 
afterglows of g amma-ray bursts have also been conducted 
with RoboPol (IKing et al.l 120141 ). More complete descrip¬ 
tions of the instrument and data redu ction p ipelin e are g iven 
in Ramaprakash et al. (in prep.) and IKing et al.l (I20141 ). re¬ 
spectively. 

The data presented here are the first obtained from an 
analysis of the instrument’s entire held of view. We present 
the observational details in Section [2] In Section [3] we de¬ 
scribe the methods used for analysing sources in the entire 
held of view. We present and discuss the results of our obser¬ 
vations in Section [4] and summarize our hndings in Section 

m 



Figure 1. A field observed by RoboPol. Each star in the field 
creates a quadruplet of images (spots) on the CCD. The central 
dark region is the mask used for lowering sky noise for the target 
at the centre of the field of view and the cross-like figure is created 
by the mask-supporting legs. 

2 OBSERVATIONS 

Observations were taken during 25 nights from August to 
November 2013, totaling around 60 hours of telescope time. 
The observations covered an area of 10 deg 2 : l = [122.6°, 
126.0°], b = [24.7°, 27.9°]. The area was initially divided 
into 275 non-overlapping helds spaced 13.2 arcminutes apart 
(slightly larger than the size of the RoboPol held of view). 
Of them, 227 were observed by the end of the period. The 
number of observations of each held ranges between 2 and 6, 
with 93% of all helds having been observed at least 3 times. 
The exposures were either 120 or 180 seconds long, with 95% 
of the exposures having lasted 180 s. All observations were 
taken in the R band. 


3 ANALYSIS 

Previous studies with RoboPol concentrated on sources ei¬ 
ther exclusively within the mask, or with the addition of 
some selected sources in the held of view around the central 
target. Although the data reduction pipeline presented by 
IKing et al] (120141 ) was designed for the entire RoboPol held 
of view, its implementation in this particular project showed 
the need for some adjustments and additions. Sources out¬ 
side the mask present a number of challenges. Some are com¬ 
mon in most polarimetric studies in the optical, while others 
are due to the particular design of the instrument. A mea¬ 
surement may be adversely affected by one of the following 
sources of systematic error: 

• Large-scale optical aberrations 

• Proximity to the mask and its legs 

• Proximity to the CCD edge 

• Proximity to other sources 

• Selection of apertures for photometry 

• Dust on optical elements 

An additional sys tematic erro r has already been iden¬ 
tified and discussed bv IKing et al] (120141 ). A rotation in the 
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Figure 2. Relative Stokes parameters across the CCD. Left: mea¬ 
sured q (top) and u (bottom). Right: residual q (top) and u (bot¬ 
tom) after subtracting the fitted model. Each square panel shows 
values on the CCD binned in 100 cells. Each cell is colored ac¬ 
cording to its average value. 



polarization reference frame of the telescope with respect to 
that of the sky causes all angles to be larger by 2.31° ±0.34°. 
All polarization angle measurements presented in this paper 
have been corrected for this. 

This section outlines the analysis of observations and 
the methodology adopted to control these systematic effects. 


3.1 Significance of measurements and debiasing 

The measurement of the fractional linear polarization ( p ) 
at the low polarization regime which is relevant for in¬ 
terstellar polarization, being always positive, is biased to- 
wards values larger than t he true (intrinsic) polarization 
(|Simmons &; Stewart 1119851 ). Thus, p measurements should 
be debiased to find their most probable intrinsic value. In the 
analysis, we consider only sources with signal-to-noise ra¬ 
tios (p/up > 2.5) so that errors are approximately normally 
distributed. The ma ximum -likeliho od es timator of the true 
value of p found by IVaillancourt I (l2006tl for measurements 
with p/a p > 3 significance is: 

Pd = vV - °p (!) 

We extend this formula to p/a p > 2.5 and use it to debias 
all measurements of p. 


3.2 Large-scale optical aberrations 

Large-scale aberrations caused by the optical sys tem are cor¬ 
rected by the instrument model, as presented bv lKing et al.l 
(I20l4l . The model is created by placing an unpolarized stan¬ 
dard star at many positions across the field of view and 
finding the best-fitting parameters that cancel the global, 
instrumentally-induced polarization and vignetting. 

The instrument model has been found to perform 



Corrected q, u across the CCD 


Figure 3. Distributions of the residuals of q and u across the 
CCD, after the subtraction of the model fit (q: dotted, u: solid). 
The vertical lines show the standard deviation of each distribu¬ 
tion. 


equally well, regardless of telescope pointing position (which 
may result in different telescope stresses) and after multiple 
removals and re-installations of the instrument on the tele¬ 
scope. The set of models that were created for these tests 
have been combined into one wi th improved perfo rmance 
with respect to that presented bv lKing et al.l (12014 '). Below 
we estimate the systematic uncertainty that remains after 
the model correction. 


3.2.1 Systematic uncertainty from model residuals 


Fig. [5] shows the uncorrected (left) and corrected (right) q 
(top) and u (bottom) values across the CCD derived by 
using this combined model. The data are binned in 100 cells 
of 39 arcsec width and the mean value is plotted in each one. 
On average, 2.4 star measurements contribute to each cell. 
The residuals appear to be homogeneous across the CCD. 

The distributions of residual q and u of the combined 
model are shown in Fig. [jj] Vertical lines show the stan¬ 
dard deviation of each distribution (a qtTeB = 0.0034, <j u , re8 
= 0.0031). Statistical errors of measurements of unpolarized 
standards are an order of magnitude lower than these stan¬ 
dard deviations, thus their contribution to this scatter can 
be ignored. Therefore, we take the systematic uncertainties 

in q and U to be Cq^sys — Uq. res, &u,sys — &V ,. res • 

From now on, in order to estimate total uncertainties 
in q and u, we add statistical and systematic uncertainties 
in quadrature, 


a 


2 

Q 


_ 2 
Uq,sys 


2 _ 2 
O'u O'u, sys 


i 2 

' ®q, stat 

(2) 

i 2 

i ^-LijStat* 

(3) 


The total uncertainty in q and u can be propagated to 
find the total uncertainty in fractional linear polarization 
(p) and electric vector position angles (EVPA or using 
the equations: 


<=^/q 2 + u 2 


q 2 a 2 + u 2 al 
q 2 + u 2 


(4) 
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Figure 4. EVPA versus fractional linear polarization of stan¬ 
dard stars. Left: BD+59.389, right: VI Cygl2. Literature values 
are shown by stars (see references in Tabled and circles are mea¬ 
surements outside the mask. Error bars include the statistical and 
systematic uncertainties added in quadrature. 
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Figure 5. Probability distribution of the intrinsic value ofp of un¬ 
polarized standards, given the measurement in the literature (red) 
and our own (black). Left: HD151406,right: HD212311. There are 
two black lines (measurements) of HD212311. 
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Assuming low polarization the expression for <7 XjSys can 
be written as: 
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surements are consistent with the literature within the errors 
(which include both the statistical and systematic uncer¬ 
tainties discussed above). Measurements ofp have not been 
debiased. 

In the case of unpolarized stars, biasing of p measure¬ 
ments is very pronounced and the interpretation of the mea¬ 
surement uncertainty is not straightforward. To facilitate 
comparison of our measurements with literature values, we 
plot, in Fig. 0 the probability distribution (likelihood) of 
the intrinsic (true) value of p, given the literature measure- 
ment (red) and our own (black). This likelihood function (see 
lYaillancourt I (120061 ). equation 8) takes into account that the 
measured values of p follow a Rice, rather than a normal, dis¬ 
tribution. In calculating the likelihood function we have used 
a total uncertainty obtained by adding statistical and sys¬ 
tematic uncertainties in quadrature, as in equations © and 
©. In both cases, our measurements are consistent within 
uncertainties with the literature measurements. There are 
two measurements (black lines) of the standard HD212311. 
For unpolarized standards the EVPA does not carry mean¬ 
ingful information, as can be seen by substituting a p /p > 1 
into equation ©: > 30°. 


3.3 Proximity to the mask, legs and CCD edge 

The mask and its supporting legs cast shadows on specific 
regions of the CCD rendering them unusable. Therefore, 
sources that happen to fall in the shadow of the mask legs 
or within 155 pixels (radially) of the mask centre are not 
considered in the analysis. 

Sources falling very close to any of the CCD edges are 
very likely to suffer partial photon losses. Also, light reaching 
these areas is subject to large optical distortions. Since the 
typical separation between a pair of the four images is 100 
pixels, we reject any spot within 100 pixels of the edges from 
the analysis. 


3.4 Proximity to other sources 


Sometimes images from different stars happen to fall within 
a few pixels of each other on the CCD. Since the typical 
diameter of a spot is 8 pixels (3.2 arcsec), photons from both 
spots are blended, as shown in Fig. [6] The relative Stokes 
parameters are computed using the following equations: 


Ni - No /4(N?a 2 + N 2 a J) 

q m+No’^ V (Ao + 2Vi) 4 


(7) 


3.2.2 Performance: Standard Stars 

To assess the accuracy of the instrument model, measure¬ 
ments of stars with known polarization values were taken 
and were then compared to the literature. During the two 
observing seasons, a number of standard stars (different 
from the ones used for the model calculation) were observed 
throughout the field of view. Catalog measurements as well 
as the band in which they were taken are shown in Table [T] 
Measurements of the unpolarized stars in the R band could 
not be found, so those in other bands are quoted. 

Fig. [1] presents RoboPol measurements of polarized 
standards (denoted by circles) and their literature values 
(stars) on the EVPA - polarization fraction plane. All p mea- 


N 3 - N 2 1 4(V|o-2 + Afcrf) 

N 3 + N 2 ,tTu \ (N 2 +N 3 ) 4 


(8) 


where Ni is the number of photons in the i th spot and m is 
the uncertainty that results from the photon noise. There¬ 
fore, overlapping of spots causes an artificially large differ¬ 
ence in intensity of one pair of spots belonging to each af¬ 
fected star. The two point spread functions (PSFs) cannot 
be de-blended, since the pipeline performs aperture photom¬ 
etry. Typically, this contamination results in erroneously 
large degrees of polarization (but not necessarily, this can 
vary based on the relative brightness of the sources involved) 
and, most notably, regular EVPAs (0°, ±45°, ±90°). This 
follows from the definition of the EVPA, equation ©. If one 
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Table 1. Polarization standard stars shown in Fig. [4] 



BD+59.389 

VICygl2 

HD151406 

HD212311 

p (%) 

6.430 ±0.022 
6.43 ±0.13 

7.893 ±0.037 
7.18 ±0.04 

0.085 ±0.041 

0.034 ±0.021 
0.02±0.021 

0.045 

x(°) 

98.14° ±0.10 
96.0 ° ±0.6 

116.23° ±0.14 
117° ±1 

-2° 

50.99° 

36.2° ±51.3° 
10.4° 

Band 

R 

R 

no filter 

V 

Reference 

1 ,5 

1, 3 

2 

1,4,5 


1 - ISchmidt, Elston &; Lupie I (1 19921) , 2 - iBerdvugin, Piir ola &; Teerikor pTl (|2014| ) , 3 - iBailev Sz Houghl \ 19821) , 4- iHeiles I d200Ql ) , 5 - 

lEswaraiah et al. I (1201 li) 



Figure 6. An example of sources that are affecting each other’s 
measurements due to their proximity (circled spots). The posi¬ 
tions of the stars (centres of quadruplets) are shown with crosses. 


Figure 8. Distributions of ratios of the FWHM between vertical 
(black) and horizontal (red) spots for a number of fields. Vertical 
lines mark the area that contains acceptable ratios. 


of the vertical images of the star is artificially brightened, for 
example N\ N 2 , then \u\ \q\ => X _■ y ±45°. Whereas if 

one of the horizontal images is affected by a nearby source, 
e.g. N 2 <C N 3 , then |q| <C |u| ^ 0°,±90°. Fig [7] (left) 

shows the EVPA versus fractional linear polarization of all 
sources with at least 2 measurements found in the Polaris 
Flare field (5172 in total). Measurements of p > 20% are 
clustered around regular EVPAs - the clear signature of 
nearby star contamination. 

We remove such sources from the analysis in the follow¬ 
ing way. If any of the four spots of a star suffers from con¬ 
fusion with another spot then we flag it as nearby contami¬ 
nated. This flag applies if a source exists within 3 x FWHM 
of a star spot. In cases where the spots of two stars happen 
to fall exactly on each other and are identified as a single 
source we check if any spot is assigned to more than one 
star. The effect of removing contaminated sources from the 
final catalog can be seen in Fig. [3 (right). All but two mea¬ 
surements with p > 20% were caused by proximity to other 
sources. 

Stars that are affected by reflections, and even other 
close-by stars in the case that the previous check fails, can 
be removed by checking the ratio of the FWHM between 
two pairs of spots. The distribution of these ratios for all 
stars found in all observed frames is shown in Fig. [8] We 


discard measurements lying outside the range 0.87 to 1.15 
(vertical lines). 

3.5 Aperture optimization 

The RoboPol pipeline performs aperture photometry to 
measure the intensity (photon counts - N) of each spot. It 
then uses these values to calculate the Stokes parameters 
as shown in equations 0 and ©• Photometry measure- 
men ts are greatly affected by the choice of aperture size 
fe.g. lHowell II 19891 1. If the aperture is too large the value ob¬ 
tained suffers from background contamination and the sig¬ 
nal to noise ratio is decreased. On the other hand, if the 
aperture is too small only a fraction of the total flux is mea¬ 
sured. This is not a problem if the same fraction of photons 
are counted, since polarimetric measurements depend on the 
relative brightness of two spots. If, however, the PSFs of two 
spots belonging to a source are different, then the fraction of 
the total flux measured is not the same and this introduces 
artificial polarization. 

A number of circumstances may affect the PSF of the 
four spots of a source. Bad seeing or weather conditions 
(wind) sometimes cause sources to appear elongated instead 
of round. Also, the optical system of the instrument may 
distort the shape of the PSF and mainly the wings of the 
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Figure 7. Fractional linear polarization versus EVPA for stars in the Polaris Flare field. Left: Measurements at regular angles (0°, ±45°, 
±90°) are caused by nearby contamination as seen in Fig. [6] Right: Measurements that survive after the removal of stars that suffered 
this contamination. Most remaining measurements of p > 5% are caused by other systematics. 


profile. Typically, bright stars (whose wings are more promi¬ 
nent) are affected more severely than faint ones. Therefore, 
it is essential that photometry be performed with an aper¬ 
ture optimized for each source. Also, the complexity of the 
optical system introduces some asymmetries in the PSFs of 
the vertical and horizontal images of a star. Consequently, 
photometry must also be optimized for each of the four im¬ 
ages of a star. 

We created a simple aperture optimization a lgorithm as 
an ad dition to the original pipeline, presented in [King et all 
d2014l) . Photon counts are measured within a circular aper¬ 
ture centred on each spot, while the background level is es¬ 
timated within an outer concentric annulus that is sepa¬ 
rated from the inner aperture by a gap. The diameter of the 
background annulus is a constant multiple of the aperture 
size. The constant is different for faint and bright sources 
as the latter have more extended wings, so that the annulus 
does not contain any light from the source while retaining 
the smallest possible distance from the source for the back¬ 
ground estimation. 

By measuring the background subtracted photon flux 
within increasing apertures we create a growth curve for 
each spot. Each of the four growth curves of the source are 
fitted with a fourth degree polynomial, P(x) (no errors are 
accounted for in the fit). The size of the aperture at which 
the normalized photon flux saturates is the optimum. To 
locate it in practice we look for the aperture size at which 
the rate of photon flux increase has reached some small value 
A. In other words, the optimal aperture is the root of the 
equation: 

f = AP (*)' < 9 > 

An example growth curve of one of the images of a star 
is shown in Fig. [3] (circles) along with its polynomial fit (solid 
line). The dashed vertical line shows the optimal aperture 
found by solving equation ©■ This aperture is used to mea¬ 
sure the photon counts and noise ( N , a) of this spot. The 
optimization is used for all data, including those collected 
for the instrument model calculation. 



Figure 9. Growth curve of one of the images of a star (circles 
show the number of background subtracted counts for each aper¬ 
ture size). A fourth degree polynomial is fit to the points. The 
optimal aperture is shown with the dashed line. 

The choice of the value for A 

To calibrate equation ([9]) and determine the best value of A, 
we created growth curves of polarization standard stars that 
were routinely observed in the field and measured their frac¬ 
tional linear polarizations and angles using all the different 
aperture sizes. Fig. 1101 shows the fractional linear polariza¬ 
tion (top) and EVPA (middle) measured for VI Cygl2 with 
different apertures. As the aperture increases, these quanti¬ 
ties saturate at some value consistent with those found in 
the literature (gray bands). As aperture size continues to 
increase, the signal-to-noise ratio worsens and also nearby 
sources may affect the measurement. The point on the hori¬ 
zontal axis after which saturation occurs is the optimal aper¬ 
ture for this star. The parameter A was selected so that it 
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reflects this transition. The bottom panel of Fig. 1101 shows 
the four growth curves of VI Cygl2 and the corresponding 
polynomial fits. The vertical line shows the aperture that 
was chosen as optimal. 

Because the standards observed with RoboPol are 
bright (typically 9-11 mag) we needed to extend this sample 
to stars of lower brightness. We selected 6 stars that were al¬ 
ready observed in the field and observed them in the mask. 
We used the values found in the mask to define the opti¬ 
mal aperture for these sources when observed in the field. 
Finally, we optimized the parameter A so that it yields an 
accurate optimal aperture for most of the stars (both these 
6 as well as the standards): A = 0.02. 

3.6 Detection of dust specks using flat field 
images 

The design of RoboPol does not allow us to correct science 
images for irregularities in transmission and uneven sensi¬ 
tivity throughout the field in conventional ways (e.g. by di¬ 
viding pixel-by-pixel by a flat field image). Because both 
sets of orthogonally polarized beams are projected on the 
same CCD, when recording an extended image (such as a 
flat) each point on the CCD is exposed to four rays tracing 
four different optical paths through the instrument. In con¬ 
trast, the photon counts we would like to correct (i.e. each 
of the four images corresponding to a point source) arrive 
on the CCD through a single optical path since each beam 
corresponds to a different orientation of the plane of polar¬ 
ization. Light from the sky against which they are projected 
still arrives through four paths for each pixel, but at differ¬ 
ent ratios, since the polarization of the sky differs between 
the moments of science and flat field image acquisition. This 
makes ordinary flat fielding impossible. 

The global non-uniformity of the field (caused by vi¬ 
gnetting) is corrected by the instrument model as described 
in Section [322] Small-scale non-uniformities cannot be cor¬ 
rected for, but they can be identified on flat-field images. 
Stars that happen to be affected by these small scale vari¬ 
ations must be excluded from the analysis. An example is 
shown in Fig. If II where the crescent pattern produced by a 
dust speck is clearly visible in the exposure and coincides 
with one of the four spots of a star (circled in white). 

We process flat-field images obtained in the evening 
and/or the morning of each night in the following way: we 
create a master flat by normalizing separate shots and tak¬ 
ing the median. After that we fit a third degree polynomial 
to the master flat and subtract its value from each pixel, 
thus removing the large scale vignetting in the flat image. 

At the position of each spot in the science image, we 
calculate the mode value (Fbgr) and standard deviation of 
counts (<7b g r) of pixels on the flat-field image that fall within 
an aperture with diameter equal to the background annulus. 
In principle, by comparing these quantities on all four spots 
of a star we can determine whether any of them has fallen 
on a dust speck, since this would cause significant variations 
in Fbgr and cj i,gr. 

To establish a set of reliable criteria that can identify 
most, if not all, dust-contaminate d stars we analysed data of 
the Be X-ray binary CepX4 (e.g. llJlmer et al. Ill973f ) , which 
is one of the mos t crowded fields observed with RoboPol 
(jReig et al. 1120141 ). We constructed a number of different 




Figure 10. Top: Fractional linear polarization of the standard 
star VI Cygl2 measured with different aperture sizes (multiples 
of FWHM). The horizontal gray band shows the literature value 
=blcr. The position of the vertical gray band shows the mean of 
the four optimal apertures while its width represents their scatter. 
Middle: EVPA measured with different aperture sizes. Bottom: 
Background subtracted number of counts with different aperture 
sizes (growth curve). The growth curves for all four of images of 
the star are shown, along with a fourth degree polynomial fit. 



Figure 11. Crescent dust patterns on the CCD. One of the four 
images of a star falls on a dust pattern (in the white circle). 
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quantities with the information from the flat field image. 
Those that proved most useful in revealing the effect of dust 
contamination were the following: 

• difference between the irbgr of a star’s vertical (horizon¬ 
tal) Spots (Ac 7bgr,v , AtTbgr,h) , 

• difference between the background value of a star’s ver¬ 
tical (horizontal) spots (AFb g r,v, AFbgr.h), 

• maximum <Tb gr (among four spots), 

• minimum <Tbgr (among four spots). 

The distributions of all 6 quantities are shown in Fig. 1121 
These quantities are measured in units of normalized counts 
in the processed master flat image. The outliers of these dis¬ 
tributions are stars that coincide with the most prominent 
dust specks. According to these distributions we selected the 
thresholds depicted by vertical lines in Fig. 1121 

Using these criteria we manage to eliminate only stars 
that are affected by the most obvious dust shadows. A 
more sophisticated analysis is needed to identify more subtle 
anomalies on the CCD. 


3.7 Statistical assessment 

The standard observing strategy in the RoboPol project is 
to obtain multiple exposures of the same field. We are thus 
able to use the stability of the measurements (in a statis¬ 
tical sense) to further reject stars with unreliable polariza¬ 
tion measurements. One reason for turning to a statistical 
treatment of the data is that even after the first stage of rig¬ 
orous cuts described in this section, some systematic errors 
are still present. These include faint dust specks, reflections 
from bright stars, and in general, sources with properties 
around the various thresholds that were used. 

First, we choose to work with stars that have reliable 
measurements of the weighted mean of p (p/crp > 2.5). The 
weighted mean is calculated by substituting into equation 
© the weighted mean q and u values of a star. 

One way to quantify the statistical significance of the 
differences between the n measurements of a star is by com¬ 
puting the reduced X 2 ( X 2 ed ) of all of its q and u measure¬ 
ments: 


n r> 

V 2 _ 1 te - <?) 

(a qJ p 

j=i 


( 10 ) 


and similarly for X 2 ed u . By placing a threshold in the value 
X 2 ed we can eliminate stars that deviate from the expected 
normal behavior. The threshold was selected so as to remove 
the tail of the distribution of all X 2 ed values of stars in the 
Polaris Flare region. The distributions of these values for 
q and u measurements are shown in Fig. [13] as well as the 
selected threshold (vertical line). 

Stars that still remain after these cuts and show signs 
of some type of contamination visible by eye on the raw 
science images were removed by hand. These include types 
of contamination already presented in this section as well 
as projected double stars, for which the analysis does not 
account. 



Figure 12. Distributions of the values extracted from the flat 
fields and used for identifying dust specks. Values on the horizon¬ 
tal axis are measured in normalized counts. 
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Figure 13. Distributions of A? , of q (black) and u (red) of stars 
with p/&p > 2.5. The vertical line shows the selected threshold. 
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Figure 14. Distribution of debiased fractional linear polariza¬ 
tions of all 648 sources resulting from the analysis. 



0 _I_I_I_I_I_I_I_ 

0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 2.0 

Ay 

Figure 15. Debiased polarization percentage vs. visual extinc¬ 
tion, Ay for all our reliably measured stars. The dashed line shows 
the maximum observable p at all Ay (p = 0.03Ay). Stars above 
the black line are marked with empty circles. 

4 RESULTS AND DISCUSSION 

The analysis provides us with 648 stars with reliable p and 
X measurements. They are presented in the online table ac¬ 
companying this paper (Table [2]). 

The distribution of debiased fractional linear polariza¬ 
tions of all these sources is shown in Fig. [14] The me¬ 
dian of the distribution is at 1.3%. Fig. QO shows the 
debiased polarization percentage a gainst visual extinction 
derive d by the SDSS colors as in ISchlaflv fe Finkbeiner I 
(1201 ill . The dashed line shows the empirically determined 
upper limit in polarization at a given Ay. p/Ay = 0.03 
dSerkowski, Mathewson fe Ford I 1 197.51 1. We mark sources 
above this limit with empty circles and use this line as a 
threshold. Sources above the line are considered separately 
as their polarizations may have an intrinsic contribution. 

In order to construct the polarization map of the re¬ 
gion we transform all EVPAs (measured with respect to the 
North-South celestial pole direction) into galactic angles ac¬ 


cording to lStephens et al~l (]201llfl . We plot the polarization 
segments of all stars below the pd — Ay line of Fig. [15] at 
each star p osition on the Herschel 250 /im image of the Po¬ 
laris Flare dMiville-Deschenes et al.~ll2010l l in Fig. 1161 The 
length of each segment is proportional to the debiased p 
of the star, calculated using equation ©■ The most strik¬ 
ing feature of the polarization map is the extended ordered 
pattern at large longitudes. In this region the plane-of-the- 
sky magnetic field appears to be oriented in approximately 
the same direction as that of the faint striations seen in 
dust emission. A very similar pattern has been observed in 
the m uch denser Taurus Molecular Cloud dChapman et al.l 
hoidf . Segments at the largest longitudes are mostly par¬ 
allel to lines of constant longitude, following the projected 
curvature of a vertical cloud structure that is partly cut-off 
by the map edges. A border appears to exist, spanning the 
diagonal of region (124°, 125°), (26°, 27°). Segments below 
this virtual line form a loop, or eddy-like feature centred at 
(124°, 25.5°) that covers latitudes down to 24.5° and lon¬ 
gitudes down to 123°. In the south, segments that are pro¬ 
jected on the dense filamentary region, also known as the 
MCLD 123.5+24.9 cloud, appear to be parallel to the axis 
of the filament and its surrounding less dense gas. A detailed 
quantitative comparison of the magnetic field as revealed by 
the map with the dust column density from Herschel will be 
presented in a follow-up paper. 

The general structure of the plane-of-the-sky magnetic 
field in this cloud agrees qualitatively with that inferred 
from the polarized emission s een b y the Planck satellite 
(IPlanck collaboration Int. XX 1120141 ). Even though the res¬ 
olution of the presented map does not allow for a detailed 
comparison, the orientation of the ordered east part is in 
fair agreement with that seen in our map. Also, the central- 
southwest part in the Planck map does show a discontinuity 
of the projected field orientation that could be a sign of the 
loop that we observe. 

The proximity of the cloud suggests that the level of 
contamination by dust foreground to the cloud is insignifi¬ 
cant. Stars lying in front of the cloud will most likely exhibit 
very low polarization (<C 1%) and so would not comply with 
the p/a p threshold, thus they would not affect the map. 

The distribution of stars for which we have reliable po¬ 
larization measurements is not uniform. Segments at higher 
galactic latitude and longitude are denser than at the lower 
part of the map. Fig. [TT] shows the number of stars in the 
map binned across the entire observed region. The bin size 
corresponds to that of the field of view. The brighter regions 
(containing more stars per bin) are in the area with ordered 
plane-of-the-sky magnetic field. This non-uniformity is not 
due to variations in the stellar density across the observed 
region. It appears as a result of the p/a p cut. We find no cor¬ 
relation between this pattern and observing conditions (i.e. 
seeing, elevation, moon phase). For all fields with a given 
number of surviving stars (N s ), we calculate the mean ex¬ 
tinction (Ay). There is a clear correlation between the two, 
as can be seen in Fig. [TH] We find that the Pearson correla¬ 
tion coefficient between these two sets is 0.59. We therefore 
conclude that this effect is most likely not the result of some 
systematic error, but that of the cloud properties. A possi- 


1 see erratum published in MNRAS. 
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Figure 16. Polarization segments over plotted on top of the Herschel 250 (im image of the Polaris Flare. The length of each segment is 
proportional to the debiased p of the star. The horizontal segment at 1 = 121° is for scale. The blue star marks the position of the North 
Celestial Pole. 


bility that could give rise to this effect is a magnetic field 
whose direction changes from mostly on the plane-of-the- 
sky in the upper left part of the map, to having a more 
pronounced component along the line of sight towards the 
lower right. 

Potentially intrinsically polarized sources 

We plot the polarization segments of sources above the 
dashed line of Fig.[l5]separately in Fig.[l9]to be easily distin¬ 
guished from those whose polarization is primarily affected 
by the magnetic field of the cloud. The length of the polar¬ 
ization segment of each star is proportional to its debiased 
p. The orientations of some segments are correlated with the 
general direction of the plane-of-the-sky magnetic field map 
of Fig. 1161 This is not surprising since the p-Ay line is em¬ 
pirical. Therefore our choice of setting a threshold based on 
that line is conservative. 

We investigate the possibility that a number of the 39 
sources falling above the pd — Ay line in Fig. [15] could be 
quasar candidates. Multi-wavelength data in this region are 
sparse, so cross-correlations with our sample were not partic¬ 
ularly fertile. The low resolution of radio data renders direct 
identification of optical counterparts impractical. For most 
highly polarized sources (over the Pd~Ay line) we only man¬ 
aged to find data from the USNOB and 2MASS catalogs. 


iKouzuma fe Yamaoka I (l201dl i presented the color proper¬ 
ties of quasar and AGN candidates in the 2MASS catalog. 
They demonstrated that candidates can be found preferen¬ 
tially at certain regions of color-color diagrams. Only one 
of our sources seems to marginally fit into this category. It 
should be noted, though, that these values have not been 
redshift-corrected. 


5 SUMMARY 

We have presented optical linear polarization measurements 
of stars projected on the Polaris Flare field. These measure¬ 
ments reveal the plane-of-the-sky magnetic field structure of 
the cloud. The observations span about 10 deg 2 of the region 
and have been conducted with the RoboPol polarimeter in 
the R-band. We presented adjustments to the automated 
data reduction pipeline that were necessary for the analy¬ 
sis of sources in the entire 13 ; x 13 7 field of view. We have 
investigated possible sources of systematic errors and have 
presented our methods for correcting for each one. 

We have produced a map of 648 polarization segments 
showing the magnetic field structure of the cloud as pro¬ 
jected on the plane of the sky. The median debiased p 
is 1.3%. The projected field shows a complicated, ordered 
structure throughout the map. To the top left part of the 


























Optical polarization map of the Polaris Flare 11 


0 ) 

3 


u 

o 

£■ 


o 



bo 

3 

'> 

V 

H 

3 

c n 

<4^ 


o 


33 

a 


3 




Galactic longitude (°) 


Figure 17. Number of stars in the map per field across the sky. 
The size of the bins corresponds to that of the field of view. The 
non-uniformity is a result of the p/<r p cut. 
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Figure 18. Number of stars (Ay) in Fig. 1161 versus mean Ay in 
all fields with N s . 


map, the field is aligned with the striations seen in dust 
emission. The bottom right parts show the presence of an 
eddy-like feature spanning roughly 2 degrees in diameter. 
Our results compare well with the Planck map of polarized 
emission of the cloud. The distribution of stars with reliable 
polarization measurements across the field is not uniform, 
with most stars lying in the top left of the region. This is 
most likely due to the intrinsic properties of the magnetic 
field structure. 
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Table 2. Reliable polarization measurements in the Polaris Flare region (full table online). 


R.A. 

Dec 

K°) 

b(°) 

Pd (%) 

<T p (%) 

x(°) 

°x(°) 

<W°) 

53.78514 

87.71787 

124.58667 

25.39440 

1.4 

0.5 

44 

7 

4 

67.42559 

88.23095 

124.53694 

26.09361 

0.8 

0.3 

36 

8 

162 

66.70140 

88.02690 

124.70368 

25.95310 

1.0 

0.3 

-14 

8 

112 

75.26493 

88.39611 

124.51862 

26.37650 

2.3 

0.6 

65 

7 
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